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Abstract. Results are presented from the first detailed spectroscopic observations of the recently identified inter- 
mediate polar RXJ0028.8+5917/V709 Cas. The study of the emission line radial velocities allows us to remove 
the uncertainties on the different aliases of the orbital period and a best value is found at (0.2225±0.0002) day. It 
is also found that the system shows significant EW~ (2-4) A broad absorptions affecting the Balmer lines from 
Hs to FL3. These broad absorptions are interpreted as the contribution of an underlying white dwarf atmosphere. 
The characteristics of the absorptions are found to be consistent with a DA log g = 8 white dwarf at a temperature 
Oh' of ~23 000K, contributing ~17% (at 4500 A) to the overall flux. This is the first direct detection of a white dwarf 

in an intermediate polar system. The absence of significant Zeeman splitting indicates a magnetic field lower 
j.^ ■ than 10 MG, confirming that, at least in some cases, intermediate polars have weaker fields than polars. Different 

' possibilities are discussed to explain the substantial contribution of the white dwarf to the overall flux. 
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1. Introduction teracting binary systems, their emission from the X-rays 

to the infrared is generally totally dominated by the re- 

Cataclysmic Variables (CVs) consist of a white dwarf ac- leage of thg gravitational energy of the accreted matter in 

creting matter from a low mass companion filling its Roche j- i j / j.- ^ j ± u j.i 

° r ° a disk and/ or an accretion column and not by the cmis- 

lobe which, in some cases, can be an evolved secondary. gion of ^ sb ^ gx components (white dwarf and low _ mass 

Among them, polars and intermediate polars, or magnetic companion) In polarS; some hints of the stellar character _ 

CVs, are a sub-class where the magnetic field of the white • ,• i , i, j j • • j ri 

' ° istics can however be obtained during the period ot low 

dwarf strongly affects the mass flow within the binary sys- accretion rates which is a typical feature of these systems . 

tern (see reviews by Cropper 1990 and Patterson 1994). During thege 1qw detected photospheric lines 

In polars (also called AM Her stars), the field is strong c£m bg uged t(j constrain the tempera ture and magnetic 

enough to lock the white dwarf spin in synchrony with fidd q{ the white dwarf (gee Beuermann 1996; Sion 1999) 

the orbit and to control the flow into an accretion column. , , t d Hi«ranre of the cnmnaninn Csee 

do W t,ll do 111L. Ildl III L- dlRl lllh LclllCL. L)l L11L. L-UllllJdjlllUll I bLL 

In Intermediate Polars (IP, previously DQ Her stars), the CrQpper ^ w&mer igg5a) posgibly due tQ & different 

field is thought to be weaker, leading to a desynchronisa- evolutionary stage and to a higher accret ion rate, among 

tion with a white dwarf spin faster than the orbit. The lpg nQ guch ^ stateg are knQwn tQ exigt &nd therefore 

accretion may then take place either through an accre- yery j.^ ig knQwn &bout the characteristics . Thc 

tion disc (disk- fed), or directly from the inner Lagrangian , , , nn , ^uji+j 

, , / only cases where a stellar spectrum has been detected are 

point (diskless) or via a combination of both (overflow- for ^ yery pecuUar long period AE Aqr &nd QK Per gys _ 

disk) (see review by P atterson 1994). Since CVs are in- tems where thg secondary is an evolved star (Welsh et aL 

Send offprint requests to: J.M. Bonnet-Bidaud 1993 ' Rehisch 1994). 

* Based on data collected at the Haute-Provence The X-ray source RXJ0028.8+5917 was detected in 

Observatory OHP (France) the ROSAT All Sky Survey and identified as an inter- 
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mediate polar on the basis of its spectral characteris- 
tics with a hard X-ray (~10 keV) component and a Nj? 
~ 10 21 at.cm~ 2 column density (Haberl & Motch 1995). 
Further pointed observations by the ROSAT satellite led 
to the detection of an X-ray pulsation at 312.8 s with 
a 40% amplitude and a probable identification with co- 
incident catalogued sources from previous different sur- 
veys 4U0027+59, 3A0026+593 and 1H0025+588 (Haberl 
& Motch 1995). Subsequent optical observations identi- 
fied the source with the m„ ~ 14 variable star, V709 Cas, 
and line velocity variations reveal two possible orbital peri- 
ods P = 5.4±0.2h and P = 4.45±0. 15 h (Motch ct al. 1996, 
hereafter M96). Optical oscillations at the X-ray spin pe- 
riod and sideband have also been reported (Kozhevnikov 
2001). A ROSAT observation in February 1998 showed a 
peculiar double-peaked X-ray pulsation with peaks sep- 
arated by one third of the spin period and the source is 
suspected to be fed by accretion through a disk seen at low 
inclination angle and with a relatively weakly magnetized 
white dwarf (Norton et al. 1999). A more detailed spec- 
tral analysis has been provided by recent BeppoSAX and 
RXTE observations (de Martino et al. 2001a, submitted) 
We report here on spectroscopic observations of the 
intermediate polar RXJ0028.8+5917/V709 Cas, which 
clearly shows the evidence of broad Balmer absorption 
lines, superimposed on the usual IP emission line spec- 
trum. These lines are attributed to the white dwarf at- 
mosphere and are used to constrain for the first time the 
temperature and magnetic field of the white dwarf in an 
intermediate polar system. 



2. Observations and Reduction 

Spectroscopic observations of V709 Cas were carried 
out at the 1.93m telescope of the Haute-Provence 
Observatory (OHP), equipped with the Carelec spectro- 
graph (Lemaitre et al. 1990). A 130 A/mm grating was 
used to obtain low resolution (FWHM « 7 A) spectroscopy 
over a wide (3600 - 7200 A) range. Observations of (2-3) 
hour duration were performed in three consecutive nights 
and an additional one separated by one day (see Table 1). 
To avoid interruption in the orbital coverage, wavelength 
and photometric absolute calibrations were performed at 
the beginning and end of the run. A careful check for the 
stability of the wavelength calibration during the run was 
made, using the strongest sky emission lines present in 
the spectrum, such as the O [I] line ( 5577.34 A), and small 
variations were corrected to an accuracy better than 0.1 A. 
The standard stars BD +28°4211 and Feige 110 were used 
for the photometric calibration on each night but due to 
variable atmospheric absorption, only the night of Aug. 
23 could be used to provide reliable absolute fluxes. To 
smooth out the influence of a possible modulation with 
the short (~ 5min) spin period, spectra were accumu- 
lated in either 15 min or 20 min exposures. All spectra were 
corrected for bias, flat-fielded and wavelength calibrated 
using the standard ESO/MIDAS procedures. 



Table 1. Log of observations 



Date 


Tstart(*) 


Tend(*) 


n 


T(min) 


1998 Aug. 21 


1047.47164 


1047.64034 


9 


20 


1998 Aug. 22 


1048.42757 


1048.64590 


11 


20 


1998 Aug. 23 


1049.49154 


1049.64089 


8 


20 


1998 Aug. 25 


1051.45438 


1051.64903 


15 


15 



(*) HJD (2450000+) , n number of spectra, T exposure times 



3. Analysis and results 

3.1. The two-component spectrum 

A total of 43 spectra were acquired over the 4 nights of 
observations. The spectra on the different nights are sim- 
ilar to those reported from the first optical observations 
by M96, showing strong Balmer and HeII( 4686 A) lines 
in emission, superimposed on a blue continuum. Figure 1 
shows a representative mean spectrum, averaged over one 
night. Weaker Hel and Hell emission lines arc also visible 
in the figure. The bump in the continuum, visible around 
5500 A , is due to broad variable atmospheric extinction 
due to ozone, not corrected by the applied standard ex- 
tinction curve. Absorptions around 6280 A and 6880 A are 
also atmospheric features. Absorption on the red side of 
Hel ( 5875A) is most likely due to the NaD doublet. Its ori- 
gin may be interstellar but it may also reveal a significant 
contribution from the secondary. 

The overall spectrum of V709 Cas is typical of IPs but 
a more careful inspection of the lines reveals an unusual 
feature. Beside the presence of an emission component, the 
Balmer lines show a very significant symmetrical depres- 
sion on both sides of the emission, suggestive of a broad 
absorption. This feature is absent in the HeII( 4686 A) 
line and the H Q line. The spectra obtained during the dif- 
ferent nights all show the same characteristics, though it 
is slightly more pronounced on Aug. 23. The normalized 
source spectrum is presented in Fig. 2, which shows more 
clearly that the most significant absorption is seen at H 7 
and H,5. The feature is less symmetric around H^ where 
the blue part of the absorption seems to be filled up. 

3.2. The orbital period 

To take into account the absorption features, the line char- 
acteristics (velocities, equivalent widths), have been de- 
rived by a fit with a two-gaussian function for emission and 
absorption, with the adjacent continuum being described 
by a one degree polynomial. The velocities of the strongest 
H Q emission line, measured over the 5 day interval, have 
been used to determine the orbital period. A Fourier trans- 
form (FFT) of the velocities was first used to constrain the 
period. The FFT yields two prominent peaks at periods 
P = 5.34±0.05h and P = 4.37±0.06h with smaller peaks 
at P = 3.70±0.06 h and P = 6.87±0.06 h, all one day aliases 
with respect to each other, the errors being determined by 
the width of the peak at half maximum. The two most sig- 
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Fig. 1. The mean spectrum of V709 
Cas on Aug. 23. The top curve is 
the observed spectrum showing a 
significant wide absorption affect- 
ing the Balmer lines from Ha to 
Hp . The middle curve is the resid- 
ual flux after subtracting a contribu- 
tion from a T=23 000K white dwarf 
whose flux distribution is shown by 
the bottom curve. Emission lines 
in the middle curve have been cut 
for clarity. The spectra have not 
been corrected for the telluric ab- 
sorption lines. The narrow absorp- 
tions and the broad bump around 
5500 A are atmospheric. The narrow 
absorption feature in the red side 
of HeI (5875 A) is identified as NaD 
(see text). 



Table 2. Emission line radial velocities 



Line 


7(km.s 1 ) 


K(km. S - 1 ) 


</>o 


xl 


H Q 


-64.7 (1.5) 


102.6 (2.3) 


0.000(3) 


4.0 




-56.7 (1.5) 


92.7 (2.3) 


-0.022(4) 


4.0 


H-Y 


-37.6 (1.5) 


88.9 (2.3) 


-0.024(4) 


5.8 


H a 


-51.0 (1.5) 


77.1 (2.3) 


-0.008(4) 


12.8 


HeII(4686A) 


-91.1 (1.5) 


75.5 (2.3) 


+0.044(4) 


7.2 


Errors on last digit in parentheses. 



nificant periods are in accordance with the values reported 
by M96. 

Velocities of the emission lines were then folded with 
the two most significant periods and fitted with a function 
V = 7 + K.sin[(27rt/P)-0o] with 7, K, and </>o as free pa- 
rameters, using a \ 2 minimization program. For all lines, 
the fit is significantly better at P = 5.34 h, so that the main 
alias period of P = 4.37 h can be securely rejected. For H a , 
the reduced Xr value decreases from 9.8 (P = 4.37 h) to 4.0 
(P = 5.34 h), corresponding to a significance at a level bet- 
ter than 99%, computed from a F-test. The (O-C) values 
of the sine blue-to-red crossing times, computed with the 
trial period of 5.34 h period, were then fitted to yield a best 
period value of P = (5.341±0.005) h with a corresponding 
ephemeris as T = HJD2451048.0575(2) + E*0.2225(2) d, 
where T is the blue-to-red crossing time of the emission 
lines. 



3.3. The line characteristics 

Table 2 gives the amplitude (K) and phasing (<f> ) of the 
different emission lines folded with the above ephemeris. 




4400 4600 
Wavelength ( A) 



Fig. 2. A close-up of the average normalized spectrum 
of V709 Cas showing clearly the absorptions around H5, 
H 7 and H^ (lower curve, shifted by -0.2 unit for clarity). 
Note the absence of such a feature around HeII( 4686 A) 
. The dotted line shows the contribution of the best fit 
white dwarf (see text). The upper curve is the residual 
after subtraction of the white dwarf contribution. 



All Balmer lines show velocity curves with similar phasing 
but with radial velocity amplitudes decreasing through 
the Balmer series from 103 km.s -1 (H Q ) to 77 km.s -1 
(Ha). The HeII( 4686 A) shows the lowest amplitude at 
~ 76 km.s -1 with a slight shift in phase and a different 7 
velocity. 

The velocity curves for Ha and H7 emission lines are 
shown in Figure 3. Inspection of the curves reveals a sig- 
nificant dispersion of the points obtained on the different 
nights. Though the minimum around <p ~ 0.75 is relatively 
well defined, there is a noticeable shift in velocity from one 
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Fig. 3. Radial velocity curves as 
a function of orbital phase for the 
H Q (bottom) and H 7 (middle) emis- 
sion components and the H 7 (top) 
absorption component. The differ- 
ent nights are shown by different 
symbols (Aug 21-squares, Aug 22- 
filled squares, Aug 23-crosses, Aug 
25-filled triangles) and the best si- 
nusoidal fits are shown by thin lines. 
Typical errors are (3-5) km.s -1 
(bottom), (5-10) km.s -1 (middle) 
and (100-200) km.s" 1 (top). 



0.5 



1.5 



night to another for points at the same phase, in the range 
4> ~(0.3-0.5) and (0.8-1.0). The dispersion from a pure sine 
curve is also reflected by the high values of the reduced 
Xr of the fits (Table 2). Using the sky emission lines, we 
checked that between the different nights the accuracy of 
the absolute calibration was better than 0.1 A. The shift 
in velocity from night to night is of the order of 70 km.s -1 
(1.0 A at H 7 and 1.5 A at H Q ). There is no clear change in 
the line profiles from night to night but we cannot exclude 
the existence of some asymmetry in the lines that would be 
poorly fitted by gaussian curves. Fits of the emission lines 
with an unique gaussian function yield a similar amount of 
dispersion, showing that the effect is not due to the exis- 
tence of the absorption component. The distortion is also 
present in H Q where no absorption is seen. The influence 
of the spin period can also be excluded as the exposure 
times of 20 min and 15 min were chosen to be close multi- 
ples of the spin period (3.8 and 2.9 respectively). Possible 
significances of this dispersion are considered in the dis- 
cussion. 

The velocity of the H 7 absorption component, as de- 
duced from the fits, is also plotted in Figure 3. The mea- 
sured values are however highly uncertain due to the very 
shallow and broad nature of the line and no significant ra- 
dial velocity curve can be derived. An attempt was made 
to measure the velocity changes with a cross-correlation 
method, using one of the best S/N spectra as a template. 
Though somewhat reducing the dispersion of the points, 
this method also does not reveal any regular variations. 

Table 3 gives the measured equivalent widths (EW) 
and intrinsic widths (FW at half-maximum) of the emis- 
sion and absorption components. For the Balmcr lines 
where an absorption is measured, the absorption is al- 
ways very broad (FWHM - 2500-4000 km.s -1 ) and with 
an EW nearly comparable to emission, with the exception 
of H/3 for which the EW is significantly smaller, due to 



Table 3. Line full widths at half maximum (FW) and 
equivalent widths (EW) 



Line Emission Absorption 
EWe(A) FW e (A) EWq(A) FW a (A) 



Ha 


18.0(4) 


14(1) 








6.6(3) 


12 (1) 


1.5 (0.2) 


36 (5) 




4.7(3) 


11 (1) 


3.5 (0.1) 


58 (7) 


H„ 


3.8(3) 


11 (1) 


3.0 (0.1) 


37 (5) 


HeII(4686A) 


6.7(3) 


13 (1) 







Errors on last digit in parentheses. 



the partial filling of the blue part of the absorption. The 
intrinsic widths of the absorption component (^40-60 A) 
are similar to what is observed in a spectrum of a DA 
white dwarf, with Balmer lines in absorption. 

4. Discussion 

The optical spectrum of V709 Cas , with strong (EW~4- 
18 A) Balmer and HeII( 4686 A) emission lines, is typ- 
ical of IPs. However, it also shows evidence of broad 
(FWHM~40-60 A) absorptions which were not yet noticed 
in this source nor in any other IP. Balmcr lines in absorp- 
tion are sometimes seen among classical CVs, mainly in 
nova-like systems in high states and during dwarf nova 
eruptions, which suggests that they are formed in an opti- 
cally thick accretion disk with a high mass transfer (see for 
instance La Dous 1994, Hessman 1986). Absorption lines 
have also been observed among dwarf novae in a very few 
systems where the accretion rate is low, as in GW Lib 
(Szkody et al. 2000) or with a high inclination, as in OY 
Car (Hessman et al. 1989). The other known cases are 
among polars during low states where Balmer absorption 
lines are clearly seen, significantly split by the Zeeman ef- 
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feet. They are attributed in this case to the atmosphere of 
a strongly magnetized white dwarf (see Cropper 1990). In 
the case of V709 Cas, the stability of the optical flux ex- 
cludes a dwarf nova event and points toward a white dwarf 
rather than a disk origin. The absence of He absorptions 
also does not favour a disk origin. The lack of any signif- 
icant absorption around H Q is indirect evidence that the 
spectral component responsible for the absorption has a 
rather steep continuum, contributing significantly only in 
the blue and becoming negligible in the red part of the 
spectrum by a dilution effect with the rest of the emis- 
sion. This is therefore indicative of a hot white dwarf as 
commonly seen among polars. 

We therefore attempted to describe the underlying absorp- 
tion spectrum, using a grid of WD models with LTE at- 
mospheres (Finley et al. 1997, Koester 2000, private com- 
munication). As the contribution of the white dwarf to 
the overall flux is unknown, we used an indirect method. 
The temperature-gravity domain was first constrained by 
comparing the observed intrinsic line widths with the ones 
predicted from theoretical models. The contribution of the 
white dwarf to the overall flux was then deduced from the 
measured EW of the absorption lines. 

For the determination of the observed parameters, we 
use primarily the H 7 measurements as the H«5 line is 
weaker and the Hg line seems to be contaminated by side 
emissions. Figure 4 shows the theoretical FWHM for a 
grid of LTE models of DA white dwarfs with a tempera- 
ture ranging from T = 15 000 to 45 000 K and a gravity 
from log g = 7.5 to 8.5. Assuming a gravity of log g — 8, 
the measured FWHM given in Table 3 is consistent with a 
temperature in the range T = 18 000-30 000 K with a best 
value at 23 000 K. If one allows for a full range of gravities, 
the T-log g domain is not unambiguously determined and 
the temperature is only weakly constrained in the range 
15 000-35 000 K. 

A dilution factor was defined as the ratio of the 
white dwarf emission to the global continuum flux, a= 
F w d/Fcont- For a continuum flux close to an emission line, 
a is also the ratio EW a /, s /EW t /,, where EW„(, S is the 
observed EW of the absorption (including dilution) and 
EW(i the theoretical EW value. The dilution factor and 
hence the contribution of the white dwarf to the global 
flux can then be evaluated from the line EW. Assuming 
log 5 = 8, the measured EW of the H 7 line in Table 3 yields 
a dilution factor of (0.17±0.04) at -4500 A for the range 
of temperatures determined above. Given the steep energy 
distribution of the white dwarf, its contribution is signif- 
icantly lower around H a with a —0.06, corresponding to 
an expected absorption component of EW a (, s — 1 A. This 
is negligible compared to the emission component EW em 
— 18 A, therefore explaining the absence of noticeable ab- 
sorption at H a . 

The comparison of the observed WD flux with the pre- 
dicted theoretical flux yields then an additional constraint 
on the R w d/D ratio shown in Figure 4, where is the 
white dwarf radius in units of 10 9 cm and D is the system 
distance in units of 100 pc. For log g — 8.0, R w d is —0.87 
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Fig. 4. Upper panel : The theoretical H 7 FWHM val- 
ues predicted from white dwarf atmospheric models for a 
range of temperatures and different log g (log g=7.5 up- 
per dotted, 8.0 middle thick, 8.5 lower dotted). The range 
of values consistent with the measured FWHM absorp- 
tions (51-65 A) and log g=8.0 are shown by vertical lines. 
Lower panel : the R^/D values, in units of 10 9 cm for 
K w d and 100 pc for D, as derived from a comparison of 
the observed and theoretical flux (see text) for the same 
parameter range. 



10 9 cm for a standard M/R relation (Nauenberg 1972) and, 
from the range of temperatures (T = 18 000-30 000 K) con- 
sistent with the observed FWHM, the system distance is 
then derived as D = (210-250) pc. 

The assumption of an underlying WD spectrum re- 
sponsible for the observed broad absorption lines leads 
therefore to a T— 23 000 K white dwarf contributing —17% 
(around 4500 A) to the overall flux for a system at — 230pc. 
These values are quite consistent with what can be ex- 
pected from a typical IP system. The absence of absorp- 
tion associated with the He lines suggests a pure DA white 
dwarf, at the difference of a disk emission where absorp- 
tion helium lines could be produced. Another interesting 
feature is the absence of Zeeman splitting, such as is ob- 
served in polars during low states. A significant Zeeman 
splitting is expected and observed in isolated white dwarfs 
for magnetic fields greater than 6-10 MG (see Achilleos et 
al. 1991) and is commonly present in polar systems with 
B fields in the range 10-30 MG. Such components are not 
detected here for H/i or H-, . With the same method as 
described in Bonnet-Bidaud et al. (2000), synthetic pro- 
files have been built using the Zeeman wavelengths and 
oscillator strengths tabulated by Kemic (1974). The ab- 
sorption features expected at the position of the Zeeman 
components from a magnetic field as low as 3 MG have 
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been compared to the observed spectra. We find in par- 
ticular that for the lowest tabulated B value of 3 MG, 
the Hp a components at 4832 A and 4896 A are clearly 
missing. The observations are therefore the first direct in- 
dication that, at least in some intermediate polars, the 
white dwarfs have lower magnetic fields than in polars. A 
low magnetic field for V709 Cas seems also suggested by 
the short (312 s) spin period of the system if one assumes 
that the system is at equilibrium with a spin period cor- 
responding to the Keplerian period at the Alfven radius 
(see Norton et al. 1999). 

The question arises as to why V709 Cas is, so far, the 
only system among IPs where the white dwarf is seen to 
contribute significantly to the overall flux so that the at- 
mospheric absorptions become visible. From the deduced 
best characteristics, there is no indication that the white 
dwarf is atypical. Its temperature (T^23 000K) is simi- 
lar, though at the upper edge of the values reported for 
magnetic CVs (Sion 1999). However, the gravity being 
not fully constrained; the contribution of the white dwarf 
could be indeed more important if the gravity is somewhat 
lower. This would imply a larger radius together with a 
higher temperature, as seen from Figure 4. For log<? = 7.5, 
the white dwarf contribution around 4500 A would be 
larger by a factor ~ 3.5 when compared to logg = 8.0. 
This would however also correspond to a rather low mass 
object (~0.34Mq) with a very peculiar high temperature 
(T - 35 000 K), the lower limit being at T > 25 000 K from 
the lower bound of the FWHM values. 

A sufficient contrast between the white dwarf and the 
overall flux could also be obtained if the contributions 
from the other regions in the system are significantly lower 
when compared to other intermediate polars. From the 
residual (total- white dwarf) flux at a distance of ~230pc, 
an absolute magnitude of M„~7.6 is derived, after correc- 
tion for a visual extinction of Ay~0.45 corresponding to 
the Nh ~ 8.7 1 20 at.cm~ 2 column density observed in the 
X-rays (de Martino et al. 2001a, submitted). If attributed 
to a disk, and using the crude (M v -M) relation given in 
Smak (1989) for face-on blackbody disks, this absolute 
magnitude would correspond, for a ~0.6Mq (logg = 8.0) 
white dwarf, to a M in the range ~ (1 — 2) 10 16 g.s -1 , 
depending on the exact extension of the disk. This is con- 
sistent with the range of values ~ (0.7—1.3) 10 16 g.s~ x de- 
rived from the recent X-ray spectral analysis (de Martino 
et al. 2001a, submitted). If, as suggested from an analysis 
of the X-ray pulsations (Norton et al. 1999), the disk is 
seen at a relatively low inclination angle, its luminosity 
would therefore not be abnormaly low. 

However, the possibility of a larger inclination angle 
cannot yet be excluded. The large dispersion seen in the 
radial velocity curves (see sect. 3.3) is strongly sugges- 
tive of the S-wave line distorsions introduced by a vari- 
able hot spot contribution and/or partial eclipses of the 
line regions found in high inclination systems (see Warner 
1995b). There is also some evidence from recent photo- 
metric observations of a possible partial disk eclipse (de 



Martino et al. 2001b, in preparation). We therefore can- 
not exclude that the important contribution of the white 
dwarf to the overall flux could be due to a combination of 
a low gravity white dwarf with a somewhat inclined disk. 

Additional observations are also clearly needed to see 
if the visibility of the white dwarf also depends on the 
overall accretion rate and luminosity of the source. 

Acknowledgments. We are very grateful to D. Koester 
and N. Dolcz for having provided us with detailed cata- 
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